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Targeted hyperthermia treatment using magnetic nanoparticles is a promising cancer therapy. However, the 
mechanisms of heat dissipation in the large alternating magnetic field used during such treatment have not 
been clarified. In this study, we numerically compared the magnetic loss in rotatable nanoparticles in 
aqueous media with that of non-rotatable nanoparticles anchored to localised structures. In the former, the 
relaxation loss in superparamagnetic nanoparticles has a secondary maximum because of slow rotation of 
the magnetic easy axis of each nanoparticle in the large field in addition to the known primary maximum 
caused by rapid Neel relaxation. Irradiation of rotatable ferromagnetic nanoparticles with a high-frequency 
axial field generates structures oriented in a longitudinal or planar direction irrespective of the free energy. 
Consequently, these dissipative structures significantly affect the conditions for maximum hysteresis loss. 
These findings shed new light on the design of targeted magnetic hyperthermia treatments. 

Tumour-targeted magnetic hyperthermia has recently attracted much attention 1 . Preferential accumulation 
of magnetic nanoparticles in tumour tissue is achieved by conjugating nanoparticles with tumour-homing 
peptides 2 or antibodies 3 . When the accumulated nanoparticles are exposed to a large alternating magnetic 
field, H=H ac sm(2nf' t), where H ac is the amplitude of the field,/is the frequency, and t is time, they begin to rotate 
because of magnetic torque. Simultaneously, the direction of the magnetic moment, //, in each nanoparticle 
reverses with a certain probability. Consequently, heat equivalent to the magnetic loss dissipates locally in the 
tumour tissue. If the properties of the irradiated field are limited (i.e., H ac 'f< constant) 1 to ensure biomedical 
safety, then nanoparticles that maximise the in vivo efficiency of heat dissipation, PhKHhc'J)> are required, where 
P H is the specific energy dissipation rate (specific loss power) per unit mass of nanoparticles. The actual rotations 
of the nanoparticles are disordered because the microviscosity of the local environment in cancer cells is not 
constant 4 ' 5 , and effective elasticity depends on the binding conditions between nanoparticles and membranes. To 
minimise the effect of irregular rotations in magnetic hyperthermia, two guiding principles have been proposed 
on the basis of simple models that consider a linear response of thermodynamic equilibrium states or magnetic 
field- driven reversals. 

The first guiding principle 6 is to use the relaxation loss in superparamagnetic iron oxide nanoparticles (SPIONs) 
with a sufficiently low energy barrier A U for reversal. If a linear response of their thermodynamic equilibrium state 
is considered at low H ac , the out-of-phase component of AC susceptibility can be expressed as follows: 

Z " = Z o-(2h/-t)/[1 + (27t/-t) 2 ], (1) 

where Xo is the initial susceptibility per unit mass of SPIONs. When reversal and rotation occur in a nanoparticle in 
parallel, the characteristic time t is given by the following equation: 

i-'^n'+^b 1 . (2) 

where t n is the Neel relaxation time for reversal, and t b is the Brownian relaxation time for rotation. 
Consequently, the heating efficiency Pu/(H ac 'f) = n[i 0 ^ r 'H ac for individual monodisperse nanoparticles has a 
single maximum at the peak frequency 2nf p = t _1 . For a sufficiently small SPION, t is determined only by t n 
because t n is much shorter than t b . In this case, it has been assumed that the conditions for maximising the 
efficiency are unaffected by uncontrolled rotation of the nanoparticles. 

However, in some experiments, dual peaks have been observed for the frequency dependence of oc P H / 
(H ac */) 7 ' 8 despite the prediction of a single peak at a 2nf p value of t" 1 (= t n _1 + t b _1 ). For this reason, size 
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distribution 7 or aggregation 8 of the nanoparticles was considered 
based on the linear response theory. In an earlier study 7 , the low- 
frequency peak observed for the susceptibility was attributed to 
Brownian relaxation of larger nanoparticles, while the high-frequency 
peak was attributed to Neel relaxation of smaller nanoparticles. In 
another study 8 , the low- and high-frequency peaks were attributed to 
individual and agglomerated nanoparticles, respectively. Thus the 
observed dual peaks have been theoretically explained by the coex- 
istence of two kinds of nanoparticles. In other words, these explana- 
tions are based on the assumption that a single kind of nanoparticle 
will produce only a single peak at t" 1 (=t n _1 + t b _1 ). However, this 
assumption has never been theoretically verified under a large AC 
magnetic field, where the linear response theory does not hold. 

The second guiding principle is to use hysteresis loss in ferromag- 
netic nanoparticles 9 . In mechanical models such as the Stoner- 
Wohlfarth model for single domain particles, ju is reversed in the 
time scale of Larmor precession (picoseconds) when A U disappears 
at the switching field H sw , because thermal fluctuations are not con- 
sidered. Such fast reversals are considered to dominate the response 
to high frequency AC magnetic field because the Brownian relaxa- 
tions of large ferromagnetic nanoparticles are generally slow com- 
pared with the oscillation of the field. In this case, the work done in 
one cycle is given by the area inside the hysteresis loop, C'M S 'H SW , 
where M s is the spontaneous magnetisation and £ is a coefficient 
related to the rectangularity of the loop. In the simple case of rect- 
angular hysteresis loops, £ is 0 for H ac < H sw and 4 for H ac > H sw . 
Consequently, the maximum efficiency, PnKH ac 'f) — C'M S '(H SW / 
#ac)/p> where p is the density of magnetite, is achieved when H ac is 
adjusted to H sw . Because H sw depends on the magnetic anisotropy 
field H K specific to each nanoparticle, it has been assumed that, in 
cases where reversal is much faster than rotation, the amount of 
hysteresis loss is unaffected by the inhomogeneous rotations of nano- 
particles in cancer cells. 

In recent experimental studies 1011 , the observed P H of immobilised 
ferromagnetic nanoparticles was lower than that of the same nano- 
particles dispersed in a fluid. Kim et al. 10 attributed the difference to 
variation in the rates of convective heat transfer. Muller et al u sug- 
gested that the orientation or agglomeration of the nanoparticles, or 
interaction effects, may be responsible for the observed difference. 
The orientation of nanoparticles is important because it is related to 
the magnetic torque intrinsic in rotatable ferromagnetic nanoparti- 
cles. However, there has been no theoretical study on magnetic field- 
driven reversals of ju in ferromagnetic nanoparticles with easy axes 
that simultaneously rotate under the magnetic torque. 

There are many reported inconsistencies between experimental 
results and predictions based on the above two guiding principles 
for optimising hyperthermia treatment. These guiding principles are 
based on the simple models established at the two limits: in zero 
magnetic field or at zero temperature. Under the conditions for 
hyperthermia (H ac ^0, where Tis temperature, the validity 

of the guiding principles has not been theoretically verified even for 
an ideal system of non-interacting monodisperse nanoparticles. 
Consequently, we attempted to simulate the thermally assisted mag- 
netic response of individual superparamagnetic/ferromagnetic iron 
oxide nanoparticles exposed to a large AC magnetic field like that 



used in hyperthermia treatment. The simulation was performed 
in the following two extreme cases: non-rotatable nanoparticles 
strongly anchored to structures resembling organelles, and rotatable 
nanoparticles in an aqueous phase mimicking cytoplasm. In the 
simulations, the thermally activated reversals of ju were calculated 
between the meta-stable directions. Simultaneously, the rotations of 
the spheroidal nanoparticles were computed in the inertialess limit 
(Brownian dynamics simulation), where the frictional torque always 
balances with magnetic torque and with Brownian torque (details are 
reported in the Methods section). The results allow examination 
of whether the relaxation loss for t n « t b and the hysteresis loss at 
H ac ~ H sw are independent of the ability of the nanoparticles to 
rotate under the conditions for hyperthermia treatment. 

Results 

The magnetic response to an AC magnetic field H ac sin(27r/"* t) at T = 
310 K was simulated for individual monodisperse spheroidal magnet- 
ite nanoparticles with non-magnetic surfactant layers in non-rotatable 
and rotatable situations (see the Methods section for details). Results 
are presented for the following representative nanoparticles: nearly 
spherical nanoparticles with an aspect ratio, /r, of 1.1 and an equatorial 
diameter, 2R M , of 18 nm, and elongated spheroidal particles with tc = 
1.4 and 2P M = 24 nm. The parameters of these nanoparticles are 
summarised in Table 1. The former nanoparticles with t n (H ac =0) 
of 20 ns can be considered as typical SPIONs, while the latter with t n 
(H ac = 0) of 2X 10 7 s (1 year) can be regarded as typical ferromagnetic 
nanoparticles in the frequency range of hyperthermia treatment 
((2tu/) _1 of approximately 1 us). Results for nanoparticles with other 
sizes and shapes are shown in the Supplementary Information. 

The magnetisation curves of the non-rotatable nearly spherical 
nanoparticles at low H ac (1 kA/m) are shown in Fig. 1. A linear 
response without hysteresis was observed at / = 100 kHz. Such 
superparamagnetic behaviour is reasonable because the estimated 
?N(H ac = 0) for the nanoparticles is 20 ns. Hysteresis appeared in 
the curves at/ = 1,000 kHz. As/increased further, the area inside the 
hysteresis loop grew. This area corresponds to the work done in one 
cycle. Therefore, PuKH&cf) a l so increased with/, and a single max- 
imum was observed at a peak frequency,/,, of 10,000 kHz (Fig. 1(b)). 
Figure 2(a) shows the H ac dependence using a contour plot of P H / 
(H ac */). As H ac increased,/, shifted towards higher frequencies. As 
indicated by the dashed line in Fig. 2(a), this shift can be explained by 
the values of T N (H ac ) calculated using the conventional Brown's 
equation as follows 12 : 

[TN(H ac )]- 1 =/o-(l-ft 2 ){(l + ft)exp[(-K d yA B r)(l+^) 2 ] 

» ( 3 ) 

+ (l-h)exp[(-K d V/hT)(l-h) 2 }} 

where h is fiH/(2 K d V), K d is the shape anisotropy constant, Vis the 
volume of the magnetic core, and k B is the Boltzmann constant. 
Therefore, the emergence of a single peak in P^KH^f) can De 
attributed to Neel relaxation loss, as expected for SPIONs. 

For the nearly spherical nanoparticles with a low H ac (1 kA/m), 
the magnetisation curves in the rotatable case were the same as those 
in the non-rotatable case (see Fig. 1). An equivalent maximum 
appeared in the /-dependence of PnliH^f) in the linear response 



Table 1 | Parameters of the simulated nanoparticles: including the aspect ratio k, equatorial diameter 2/?m, thickness of the surfactant layer 
3R, density of magnetite p, spontaneous magnetisation M s , anisotropy field Hk, viscosity 77, Neel relaxation time tn(H = 0), and Brownian 
relaxation time tb. The value of tn(H = 0) was estimated using equation (3), while the value of tb was estimated using equation (4). 

2/? M SR P M s H K t] t n (H=0) T B 

k (nm) (nm) (kg/m 3 ) (kA/m) (kA/m) (mPa s) (sec) (sec) 

Nearly spherical nanoparticle 1.1 18 4.5 5200 450 17 l.O/oo 2xl0" 8 8xl0" 6 

Elongated spheroidal 1.4 24 6.0 5200 450 57 l.O/oo 2xl0 7 2xl0" 5 

nanoparticle 
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(a) M-H curves (b) Magnetic loss 




f (kHz) 



Figure 1 | Magnetic response of nearly spherical nanoparticles (typical 
SPIONs). The open and solid symbols show the values of non-rotatable 
and rotatable nanoparticles, respectively, (a) Steady magnetisation curves 
with low H ac ( 1 kA/m) at various frequencies, (b) frequency dependence of 
the efficiency of heat dissipation with low and intermediate values of H ac . 

range (H ac = 1 kA/m) (Fig. 1(b)). This behaviour is consistent with 
the above assumption because the estimated t n (H ac = 0) of 20 ns is 
much shorter than t b = 8 jas (Table 1). The shift of this peak with 
increasing H ac is analogous to that in the non-rotatable case 
(Fig. 2(b)). However, another maximum of Pul{H ac -f) was observed 
at H ac = 8 — 16 kA/m and/ = 30 kHz in the contour plot shown in 
Fig. 2(b). A secondary maximum like this has not previously been 
theoretically predicted for individual monodisperse nanoparticles. 
Figure 3(a) shows the magnetisation curves calculated under these 
conditions. Unlike the non-rotatable case, an S-shaped hysteresis 
loop without remanence existed. At the same time, the mean 
orientation of the long (easy) axes of the nanoparticles showed 
butterfly- shaped hysteresis, as shown Fig. 3(b). Because such beha- 
viour cannot be explained using the present linear response theory, 
its origin is discussed in the next section from the viewpoint of the 
rotation of the long axis of SPION. 




H ac (KA/m) H ac (kA/m) 

Figure 2 | Efficiency of heat dissipation in the nearly spherical 
nanoparticles (typical SPIONs) that are (a) non-rotatable and (b) 
rotatable, where P H is the specific energy dissipation rate. Dashed lines 
represent the Neel relaxation time (27iT N ) -1 , dotted lines show the 
Brownian relaxation time (271Tb)" 1 , and solid lines indicate typical angular 
velocity, w H (H = H ac , \]/ = 7i/4)/2ti, of the rotation caused by magnetic 
torque. White lines show the thresholds for biomedical safety. Diamonds A 
and B on the white lines denote the conditions for maximum Pn/(H ac 'f) in 
the rotatable nanoparticles. 



(a) M-H curves (b) Mean orientation 




(kA/m) 



Figure 3 | Magnetic response of nearly spherical nanoparticles (typical 
SPIONs) with an applied AC field with H ac = 8 kA/m and f = 30 kHz. 

(a) Steady magnetisation curves, (b) mean orientation of the long (easy) 
axis of the nanoparticles, <cos0>. Orientations are indicated in the inset 
images. 

The magnetisation curves at / = 10,000 kHz for the elongated 
spheroidal nanoparticles, which are typical ferromagnetic nanopar- 
ticles, in the non-rotatable case are shown in Fig. 4(a). The curve was 
reversible at H ac = 20 kA/m, and hysteresis appeared in the curve at 
H ac = 26 kA/m. As H ac increased further, the area inside the hyster- 
esis loop grew. When H ac became larger than 32 kA/m, the expan- 
sion of the area was saturated and the shape of the magnetisation 
curve approached that predicted by the Stoner-Wohlfarth model. 
PuKHac'f) was almost zero at low H ac , then at approximately 
30 kA/m it began increasing rapidly with H ac , followed by a gradual 
decreases with increases in H ac (Fig. 4(b)). This behaviour depends 
weakly on the frequency (Fig. 5(a)). In mechanical models that do not 
consider thermal fluctuation, a hysteresis loop appeared when H ac 
was higher than H sw . Because H sw of ferromagnetic nanoparticles 
with randomly oriented easy axes ranges from H K /2 = 29 kA/m to 
H K = 57 kA/m, and is often close to 0.5 H K , that PnHfi^f) is almost 
independent of frequency at H ac > 30 kA/m in the non-rotatable 



(a) M-H curves (b) Magnetic loss 




Figure 4 | Magnetic response of elongated spheroidal nanoparticles 
(typical ferromagnetic nanoparticles) in high-frequency AC fields with 
f= 10 000 kHz and various values of H ac . (a) Steady magnetisation curves 
in the non-rotatable case; the corresponding curves for rotatable 
nanoparticles are presented in Fig. 6(a) and (b). Solid lines indicate 
Stoner-Wohlfarth model curves, (b) H ac — dependence of the efficiency of 
heat dissipation. The open and solid symbols show the values of the non- 
rotatable and rotatable nanoparticles, respectively. The arrows indicate the 
peak maxima of PuKH ac 'f) in both cases, and the broken line shows half of 
the anisotropic field, H K /2. The inset shows an enlarged view. 
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(a)non-rotatable 
30000 




(a) Major loops 



(b) Minor loops 



4 16 k641 

H ar (kA/m)\ 



16 |K 64 



A7 (kA/m)\ 



Figure 5 | Efficiency of heat dissipation of elongated spheroidal 
nanoparticles (typical ferromagnetic nanoparticles) that are (a) non- 
rotatable and (b) rotatable. Dashed lines represent the Neel relaxation 
time (27iT N ) _1 and the solid line indicates f p , which was calculated using 
equation (5). White lines show the thresholds for biomedical safety. 
Diamonds C and D on the white lines denote the conditions for maximum 
PuKH ac 'f) in the rotatable nanoparticles. 

ferromagnetic nanoparticles is consistent with the properties 
expected for the hysteresis loss. 

Figure 6 shows the magnetisation curves for rotatable elongated 
spheroidal nanoparticles at/= 10,000 kHz. Because the magnetic 
response slowly changed after the AC magnetic field was applied at 
t = 0, transient variations of the hysteresis loops are observed. The 
shape of the major hysteresis loop at H ac = 60 kA/m was initially 
consistent with that predicted by the Stoner-Wohlfarth model with 
randomly oriented easy axes. However, the remanence of the major 
loop gradually increased from 0.5 M s to M s . In other words, the 
major loop became squarer, and the area inside the loop increased 
with time. In comparison, the remanence of the minor loop at 
H ac = 26 kA/m gradually decreased and the area became smaller 
over time. As shown in Fig. 6(c) and (d), the long (easy) axes of the 
nanoparticles gradually turned when the variations of the loops pro- 
ceeded (see the next section for details). Consequently, the increased 
area of the major hysteresis loops and decreased area of the minor 
loops caused the maximum of PuKH^'f) to shift towards higher H ac 
compared with the non-rotatable case (see arrows in Fig. 4(b)). Note 
that reversals occurred every hundred nanoseconds (~l/f), while 
rotations took several microseconds (Fig. 6(c)). Thus, the assump- 
tion that the amount of hysteresis loss is unaffected by the rotation 
of nanoparticles when reversal is significantly faster is invalid 
for ferromagnetic nanoparticles in large AC magnetic fields at high 
frequencies. 

Before discussing this novel phenomenon observed at high 
frequencies, the other important variation in the contour plot of 
Pu/(H ac 'f) (Fig. 5) that occurred because of the ability of the fer- 
romagnetic nanoparticles to rotate at lower frequencies shall be 
examined. The maximum of PuKH^'f) shifted toward lower H ac 
below 100 kHz for the rotatable elongated spheroidal nanoparticles, 
while it stayed between H K /2 and H K in the non-rotatable case. 
Figure 7(a) shows the magnetisation curve obtained when H ac = 
16 kA/m and / = 30 kHz. The curve in the rotatable case had an 
obvious hysteresis loop with a large remanence in the steady state, 
but there was no hysteresis observed for the non-rotatable situation. 
Because H ac =16 kA/m is much smaller than H K /2, no reversals of ju 
occur at any orientation of the easy axis so hysteresis is not observed for 
the latter case. Figure 7(b) shows the variation of <cos0> in the 
rotatable case, where <cos0> (0 < 0 < n/2) is the mean angle between 
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(c) Mean orientation (d) Relaxation time 
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50 
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Figure 6 | Magnetic response of elongated spheroidal nanoparticles 
(typical ferromagnetic nanoparticles) in high-frequency AC fields with 
f = 10 000 kHz and various values of H ac . (a) Transient major hysteresis 
loops of the rotatable nanoparticles after application of a field with H ac = 
60 kA/m at t = 0. (b) Transient minor loops of the rotatable nanoparticles 
after application of a field with H ac = 26 kA/m at t = 0. Solid lines indicate 
Stoner-Wohlfarth model curves, (c) Relaxation of the mean orientation of 
the long axis, <cos 6>, (d) steady values of <cos 6> and relaxation time 
T rot . Dashed lines in (d) show H K /2 and H K ; and the solid line in (d) 
indicates the reciprocal, [co H (H=H ac sin^ = 0.43)] _1 , of the typical angular 
velocity of rotation caused by magnetic torque. Orientations are indicated 
in the inset images. 

the magnetic field and the long axes of the spheroidal nanoparticles. 
Note that <cos0> is synchronised with \M/M S \ = |cos^|, where \jj is 
the angle between ju and H. This fact indicates that the hysteresis in the 
rotatable case (Fig 7(a)) is mainly caused by the rotation of the easy axis 
where the direction of ju is fixed. Consequently, heat equivalent to the 
hysteresis loss dissipates even at H ac < H K /2. 

For Brownian relaxation, t b can be expressed as follows: 

TB = (3iyV H -(0.8 + 0.2fc)/(fc B r)), (4) 

where y\ is the viscosity of the surrounding medium, and V H is the 
hydrodynamic volume. In equation (4), the factional torque for 
spheroids, described in the Methods section, is considered. For the 
elongated spheroidal nanoparticle, (27TTb) _1 is calculated to be 
8 kHz. This value is too low to cause the nanoparticle to rotate at 
30 kHz. Therefore, Yoshida et at. 13 also took into account the rota- 
tion caused by magnetic torque, ju{t) XH{t). They concluded that the 
area of the hysteresis loop was maximised as follows: 



27r/ p = T- 1 [l + 0.07(^H ac //c B T) 2 ] 



(5) 



The location of the peak in PuKHac'f) below 100 kHz can be 
explained by this equation, as shown in Fig. 5(b). An expression that 
describes all of the variation in the position of the primary peak of 
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(a) M-H curves (b) Mean orientation 




H(kA/m) 



Figure 7 | Magnetic response of elongated spheroidal nanoparticles 
(typical ferromagnetic nanoparticles) in a low-frequency AC field with 

f = 30 kHz and H ac =16 kA/m. (a) Steady magnetisation curves, where 
the open and solid symbols show the values of the non-rotatable and 
rotatable nanoparticles, respectively. Solid lines are Stoner-Wohlfarth 
model curves, (b) Mean orientation of the long (easy) axis of nanoparticles, 
<cos#>. 



PHKH ac 'f) is desirable, and equation (5) can be rewritten as follows: 
2nf p ~ [T N (H ac )]- 1 + T- 1 [l + 0.07( A iH ac //c B T) 2 ] a5 . (6) 

This equation is an extended relationship of t" 1 = t n _1 + t b _1 for 
a large AC magnetic field. This expression is for the primary 
maximum; the secondary maximum is discussed later. The second 
term of equation (6) can be approximated to 0.1 /nH ac /(tiV H ) for 
fiHJk B T » (0.07)"° 5 and tc — 1. On the other hand, T N (H ac ) of 
ferromagnetic nanoparticles becomes extremely short only when 
the energy barrier disappears between H K /2 and H K . Therefore, the 
changeover between the two terms in the equation (6) generally 
occurs at H ac - H K /2 and 2tc/« 0.1/tff K /(2*/V H )=0.1 K d V/(rjV H ). 
For the elongated spheroidal particles (k = 1.4, iC d =16 kj/m 3 , and 
VVVr — 0.3) in a liquid phase with r\ = 1 mPa*s, the values of H ac 
and /are 29 kA/m and 76 kHz. Such a changeover around 100 kHz 
occurs for ferromagnetic nanoparticles of any size, as long as the 
conditions, K d , V/V H , and r\, are constant. We must keep in mind 
that, even when ferromagnetic nanoparticles are large enough for 
their Brownian relaxation to be negligible, the magnetic torque 
caused by the large AC magnetic field can easily rotate such nano- 
particles in the liquid phase at a time scale of microseconds. This 
knowledge is helpful when considering the frequency for hyperther- 
mia treatment, even if it is obtained for a simplified system. 

In summary, most of the simulated results, including signifi- 
cant variations for ferromagnetic nanoparticles exposed to a low- 
frequency AC magnetic field, can be explained using the existing 
models. The two essential exceptions are as follows: 

(a) a secondary maximum in the relaxation loss for SPIONs exposed 
to a low- frequency AC magnetic field, and 

(b) a shift of the maximum hysteresis loss caused by the ability of 
typical ferromagnetic nanoparticles exposed to a high-frequency AC 
magnetic field to rotate. 

These novel phenomena are discussed in detail in the following 
section. 

Discussion 

The two novel phenomena of rotatable nanoparticles in a large AC 
magnetic field described above cannot be explained by simple models 
that consider a linear response of thermodynamic equilibrium states 
or magnetic field-driven reversals. In this section, we shall further 



discuss these atypical responses from the viewpoint of the orientation 
of the long (easy) axis. First, we begin with the appearance of a 
secondary maximum of PuKH^'f) near H ac = 8 kA/m and / = 
30 kHz for the nearly spherical nanoparticles (typical SPIONs) 
where an S- shaped hysteresis loop without remanence was obtained 
(Fig. 3). We must recall that the variation in <cos0> showed 
butterfly- shaped hysteresis under these conditions. This behaviour 
explains the atypical magnetic response in the period/ -1 (33 jas) 
(Fig. 3(b)). Initially (at t = 0), no magnetisation exists because the 
occupation probabilities of ju in the two stable directions parallel to 
the long (easy) axis are equalised in a zero magnetic field. As H 
increases, the occupation probability in the more stabilised direction 
immediately increases because of reversals on a time scale of t n 
(< 20 ns). The reversed //in the stabilised direction is not completely 
parallel to H, \jj 0, and the magnetic torque /nHsmij/ turns the long 
(easy) axis towards the direction of the field. If we neglect Brownian 
torque A{t) (see equation (11) in the Methods section), the angular 
velocity of the rotation due to magnetic torque can be expressed as 
o> H (H(f), m) = W)siniA(0]/NV H -(0.8 + 0.2K)]. (7) 

Hence, co H (H(t), \j/(t)) increases in proportion to the field amplitude 
H = H ac sin(27i/f). For example, w H (H, i// = n/4) is 0.15 X10 6 rad/s 
when H (t = 1/4/ ~ 8 |is) is 8 kA/m. Therefore, rotation is not 
negligible in the peak period of the oscillations of H. Subsequently, 
H decreases to zero at t = 1/2/ ~ 17 |is, and the occupation prob- 
abilities are again equalised because reversal is rapid, so the magnetic 
torque disappears. Alternatively, the Brownian torque randomises 
the orientation of the long axis on a time scale of t b (= 8 |is). 
Therefore, competition between the magnetic and Brownian torques 
can cause the butterfly- shaped hysteresis of <cos0>. Because the 
equilibrium magnetisation of SPIONs with easy axes parallel to His 
higher than that of randomly oriented SPIONs 14 , the magnetisation 
curve shows hysteresis without remanence. Consequently, a second- 
ary maximum appears for the rotatable SPIONs even though t n « t b . 

Next, we investigate the influence of the ability of elongated spher- 
oidal nanoparticles to rotate under a high-frequency AC magnetic 
field on the shift of the maximum PuKH^'f)- As shown in Fig. 6(a) 
and (b), the magnetisation curves varied after an AC magnetic field 
was applied at t = 0. At the same time, transient variations also 
occurred in <cos0>, as shown in Fig. 6(c). In the case of the major 
loop at H ac = 60 kA/m, <cos0> gradually increased from 0.5 to 
0.95. In other words, the long (easy) axis became oriented towards 
the direction parallel to H. The characteristic time, T rot , was estimated 
to be 1.3 |is using the approximation of exponential decay. Note that 
the direction of ju is not completely parallel to H even for H > H K , 
even though //is already reversed for all of the nanoparticles. Because 
sin \jj is 0.43 when cos \\i is 0.9, a large magnetic torque can turn the 
long axis even if the magnetisation is almost saturated after reversals 
at H ~ H K . Indeed, T rot at H ac > H K is comparable to the reciprocal of 
typical values of Wn(H = H ac , sin \jj = 0.43) (Fig. 6(d)). Therefore, 
these transient variations can be attributed to the longitudinal ori- 
entation being adopted preferentially because of the magnetic torque. 

In the minor loops at H ac =26 kA/m, the remanence of the rotat- 
able nanoparticles decreased gradually with time (Fig. 6(b)) and 
<cos0> simultaneously decreased from 0.5 (Fig. 6(c)). The long axis 
was oriented perpendicular to H during this period, although the 
longitudinal orientation is preferred when the Zeeman energy is 
considered. It should be noted that the angle xj/ for ju in a stable 
direction more parallel to H is smaller than that in a metastable 
direction less parallel to H (Fig. 6(d)). In other words, the magnitude 
of the magnetic torque toward the longitudinal orientation in the 
former is weaker than that toward the perpendicular orientation in 
the latter. This difference makes the orientation of the long axis 
planar on average, because the stable and metastable states alternate 
every half period when the reversal of //is blocked in the minor loop. 
These arguments suggest that the slowing of the rotation for 
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0.5 H K < H ac < H K (Fig. 6(d)) can be attributed to compensation 
between two magnetic torques in the intermediate range. Briefly, in 
ferromagnetic nanoparticles in the aqueous phase, longitudinal or 
planar orientations are adopted, irrespective of the free energy, as 
dissipative structures under a high-frequency AC magnetic field. 
Consequently, PuKH ac 'f) increases gradually in major hysteresis 
loops and decreases in minor loops. These variations cause the max- 
imum of PuKHac'f) t° shift towards higher H ac . 

Finally, we return to the contour plots of Pn/(H ac 'f) in Figs. 2 and 
5, and discuss the effect of rotation on the design for maximising P H / 
(Hac-/). If a safety limit of H ac */< 4.85X10 8 Am'V 1 is applied 9 , 
then maximum values of PuKH ac 'f) fc> r rotatable SPIONs and fer- 
romagnetic nanoparticles are obtained at the conditions shown by 
diamonds A and C in Figs. 2(b) and 5(b), respectively. However, no 
heat dissipation occurs under the same conditions (A and C) if the 
rotation of these nanoparticles is blocked (Figs. 2(a) and 5(a)). If a 
highly amplified AC magnetic field H ac */of 1.74X10 10 Am'V 1 is 
allowed 15 , a maximum P H /(H ac */) of 3.0X 1(T 4 m^s^A" 1 
(5.2 MW/kg) for the rotatable ferromagnetic nanoparticles can be 
obtained (diamond D in Fig. 5(b)). However, PuKH ac 'f) halves when 
the rotation of these nanoparticles is blocked (Fig. 5(a)) because 
oriented structures are not formed. In contrast, condition B for the 
primary maximum of Pn/(H ac 'f) in the rotatable SPIONs remains 
the optimum condition when these nanoparticles cannot rotate 
(Fig. 2(a) and (b)). This is because the long (easy) axes of SPIONs 
are randomly oriented in rotatable SPIONs as Brownian torque has 
more effect than magnetic torque in a weak magnetic field. As 
demonstrated here, rotation generated by the magnetic torque 
caused by a large alternating magnetic field greatly affects the con- 
ditions for maximising heat dissipation in magnetic nanoparticles. 

In this study, we simulated the magnetic responses of superpar- 
amagnetic and ferromagnetic magnetite nanoparticles in a large 
alternating magnetic field. The results show that both the relaxation 
loss for t n « t b and the hysteresis loss at H ac ~ H sw are affected by the 
formation of dissipative structures because of the ability of nanopar- 
ticles to rotate. Consequently, the conditions for maximising heat 
dissipation depend strongly on the inhomogeneous microviscosity of 
the surrounding medium. 

Compared with the simplified model used for our simulation, 
actual magnetic nanoparticles used for targeted magnetic hyperther- 
mia treatment are not ideal. For this reason, the factors affecting 
more realistic situations need to be evaluated. First, the effects of 
crystalline and surface anisotropy energy are considered. In this case, 
the potential energy with respect to the direction of// is complicated. 
Even if multiple valleys appear in the energy surface, the easy axes are 
not parallel to H, because the orientations of the nanoparticles are 
randomised by Brownian torque in the liquid phase. For this reason, 
slow rotations inevitably occur after fast reversals because of the 
magnetic torque in an AC magnetic field. These rotations lead to 
secondary relaxation loss in SPIONs in a low frequency AC magnetic 
field and shift the hysteresis loss in ferromagnetic nanoparticles in a 
high frequency AC magnetic field. Next, the variation in the size and 
shape of actual nanomagnets must be considered. In this case, Neel 
relaxation times, t n , differ significantly because they depend expo- 
nentially on the volume of each nanoparticle and the shape aniso- 
tropy constant. In contrast, the dependence of frictional torque on 
the size and shape of nanoparticles is weak. Because the S-shaped 
hysteresis loop of SPIONs appears in the frequency range of rotation, 
the secondary loss peak becomes less diffuse compared with the 
primary relaxation loss peak. For ferromagnetic nanoparticles, the 
shift of the hysteresis loss at high frequencies should still be signifi- 
cant even if the size of nanoparticles is not uniform, because the 
anisotropy field is independent of nanoparticle size. Finally, the effect 
of dipole-dipole interactions is considered, because the density of 
nanoparticles accumulated in cancer cells might be inhomogeneous 
if they are trapped at specific sites. In such a case, chain structures of 



longitudinally aligned nanoparticles have been conventionally dis- 
cussed in a magnetic field, although their details are still controver- 
sial 16 . Our findings illuminate this conventional view, because, in 
some cases, formation of structures with a planar orientation is pre- 
dicted even for individual ferromagnetic nanoparticles. In future 
studies, we will clarify a variety of dissipative structures, which 
are different from ordinary chains, for interacting nanoparticles. 
As discussed here, knowledge of the heat dissipation in the non- 
equilibrium steady states of rotatable nanoparticles is essential for 
the design of targeted magnetic hyperthermia treatments using large 
AC magnetic fields. 

Methods 

Model for the simulation. The magnetic response to an AC magnetic field 
H ac sin(27i/"* t) was simulated for individual superparamagnetic/ferromagnetic 
magnetite nanoparticles in two extreme cases: non-rotatable nanoparticles strongly 
anchored to structures resembling organelles and rotatable nanoparticles in an 
aqueous phase resembling the cytoplasm. We considered the nanoparticles to be 
monodisperse prolate spheroids with equatorial diameters, 2R M , from 12 to 24 nm, 
and aspect ratios, k, between 1.1 and 1.4. Because these dimensions are smaller than 
the typical exchange length of magnetite, 27 nm, 17 all of the spins are parallel to one 
another in each nanoparticle. In other words, we can assume that each nanoparticle 
has a single magnetic moment ju = M S V (coherent rotation/ macro-spin 
approximations 14 ), where Vis the volume of each nanoparticle, [(4/3)k , kR m 3 ]. The 
magnitude of spontaneous magnetisation M s = |M S | was set to the value of bulk 
magnetite 6 , 450 kA/m, because the dependence of M s on particle size has not been 
well established 18 . 

The magnitude of ju is important; for example, \i of a nanoparticle with 2R M = 
24 nm and k = 1.4 is 1.5X 10 5 fi B . Hence, such nanoparticles aggregate because of the 
large dipole-dipole interactions between ju unless a sufficient non-magnetic sur- 
factant layer exists. The required thickness of this layer is approximately 6 nm for 
nanoparticles with 2R M = 24 nm 19 . Therefore, in our model, we used non-magnetic 
layers with a thickness 3R of 0.5 R M . Consequently, we can assume that the nano- 
particles are uniformly dispersed and do not aggregate. In this case, the typical 
distance between nanoparticles is n~ m , where n is the number density of nanopar- 
ticles. Because the actual mass fraction of nanoparticles accumulated in cancer cells 
does not exceed 1% (approximately 10 mg/mL), the magnitude of the interaction 
between nanoparticles with 2R M = 24 nm, n'fi 2 /k B , is estimated to be less than 
IX 10 2 K. Consequently, we did not consider minor effects caused by dipole-dipole 
interactions in our simulations. 

In such an individual nanoparticle, the potential energy with respect to the dir- 
ection of //, U(Q)> is given by the following equation: 

U(Q) = U d (Q) + U C (Q) + U S (Q) -fi H (8) 

where U d (Q), U C (Q), U s (£2) and —//'Hare the shape, crystalline, and surface 
anisotropy energies, and Zeeman energy, respectively, and Q is the solid angle of //. 
For spheroidal particles, the first term, U d (Q), can be described as KdVsin 2 ^ and has 
two minima separated by the energy barrier AU=K d V, where K d is the shape 
anisotropy constant and ^ is the angle between the long (easy) axis and //. The 
magnitudes of K d , given by (N e asy-Nh a rd)/V^s 2 > are 5 and 16 kj/m 3 for spheroidal 
particles with k of 1.1 and 1.4, respectively. N easy and N hard are the demagnetising 
factors for the long and short axes, respectively. In comparison, magnetite has cubic 
crystalline anisotropy 20 and an anisotropy constant K x = — 11 kj/m 3 . The energy 
landscape of cubic anisotropy (K^O) is gentle and Al/is (1/12)^^ « (1 kj/m 3 )- V 19 . 
Additionally, the effects of surface anisotropy are generally insignificant if the part- 
icles are larger than 10 nm 18 , although the rationale for this is still unknown. For these 
reasons, we assumed that the uniaxial shape anisotropy dominated. Consequently, 
equation (8) can be simplified as follows: 

U(<p,\j/) ~K d V sin 2 (p-fiH^ sin(2nft) cos^, (9) 

where \J/ is the angle between ju and H (see Fig. 3). 

Little is known about the local environments of magnetic nanoparticles accumu- 
lated in tumour tissue. For example, nanoparticles coated with dextran are completely 
immobilised in tissues 21 , whereas dextran nanoparticles appear to be mobile in cells 5 . 
For this reason, we simulated the magnetic response of the monodisperse spheroidal 
nanoparticles in the following two extreme cases: non-rotatable nanoparticles with 
randomly oriented easy axes, and rotatable nanoparticles in a Newtonian fluid with a 
viscosity rj of 1 mPa*s. 4 

Simulation of the reversal. The detailed trajectories of// in a magnetic field applied at 
an oblique angle 9 to the long easy axis of a spheroidal particle can be precisely 
simulated by solving the stochastic Landau-Lifshitz-Gilbert equations. However, we 
are only interested in the reversal of// once every microsecond because the frequency 
used for hyperthermia is limited. Therefore, we can use a well-known coarse-grained 
approach, or "two-level approximation" 14 , that considers thermally activated 
reversals between the meta-stable directions, ((f)i, i/^) and (<^ 2 > i/^)* via the midway 
saddle point at (<^ 3 , in Uq(<I),iJ/). The reversal probability from (</>!, i/^) to (0 2 > 
D12, is given by/ o _1 'exp[((7e(03, i/^ 3 ) — Uq^i, iJ/i))/k B T], while the backward reversal 
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probability u 2 \ is/o _1 *exp[(l/e(03, \j/ 3 ) — Uq((\) 2 , ^2))/k B T] y where f 0 is the attempt 
frequency of 10 9 s" 1 . 

In the simulation, the time evolution of the occupation probabilities, pi(6), p 2 (0) = 
1 — pi(0)> at the two stable directions of a nanoparticle tilted at 9 was computed using 
the following relationship: 

dp 1 /dt = V2ip 2 (0)-v 12 p l (9). (10) 

pi(6) was simply set to either zero or one when there was only one minimum in 
Uq((/>,\I/). This calculation was continued until the transient factors depending on the 
initial conditions disappeared. The time step At was typically 10" 4 //"s but was shorter 
when VyiAt (or v 2i At) became large compared with one. At each step, magnetisation 
was obtained as \{pi(9)ficos\J/i + p 2 {9)picos\J/ 2 )sm9d9. Test simulations were per- 
formed to check the validity of this method using the same parameters as in an earlier 
study 14 to calculate the reversals of magnetic nanoparticles in large AC magnetic 
fields. As detailed in the Supplementary Information, our results agree with the 
reported behaviour 14 . 

Simulation of reversal and rotation. The rotation of spheroidal nanoparticles was 
simulated in synchronisation with the reversal of ju. In Newtonian fluids, the frictional 
torque for rotation can be expressed as 6rjV H '(0.S + Q2k)' co(t), 22 where V H = [(4/ 
3)k'k(R m + SR) 3 ] is the hydrodynamic volume and at(t) is the angular velocity of 
rotation; de/dt = co{t) X e (£), where e (t) is the unit vector along the long axis of the 
spheroid; and ju{t)'e(t) = /icos</>. Under typical conditions, where rj = I mPa*s, V H 
~ 10 3 nm 3 , and co(t) ~ 1 X 10 5 rad/s, the inertia of the nanoparticle can be neglected 
(Brownian dynamics simulation). In this inertia-less limit 23 , the frictional torque 
balances with magnetic torque /t(t)XH(t) and Brownian torque A(t) as follows: 

6riVn-(0.8 + 0.2Kyco(t) = n(t)xH(t) + A(t) (11) 



<Ai(f)>=0, (12) 

<A i (t l )l i (t 2 )>=2k B Ti6riV li i0.8 + 0.2K)y3(t l -t 2 ), (13) 

where 8(tl— 12) is the Dirac delta function. 

At the beginning of the simulation for the rotatable nanoparticles, an assembly of 
randomly oriented nanoparticles was generated, where their number ensures an 
optimal compromise between calculation time and precision. Then, the time evolu- 
tion of the direction of ju and the orientation of the long easy axis were computed by 
the following steps, (i) Using equation (9), reversible variations of the meta-stable 
directions (^(t), caused by the latest changes in the field strength and direction 
of the easy axis were calculated, (ii) ju{f) at (0i(f)» ^i(f)) was reversed if x < u 12 At, but 
otherwise not. In this case, x e [0, 1] is a pseudorandom number generated by the 
Xorshift algorithm 24 . The backward reversal was computed in a similar manner, (iii) 
Substituting the reversed (or held) ju{f) into equation (11), co{t) was calculated, (iv) 
e(t) was finally computed using the relationship de/dt = co{t)X.e{t). This calculation 
was continued until transient factors depending on the initial conditions disappeared. 
In this simulation, At was typically 10~ 4 //s but was shorter unless Vi 2 At, or the 
changes in e{t) were sufficiently small compared to one. Magnetisation was obtained 
as ^(Wi/icos^i) at each step. Test simulations were performed to check the validity of 
this method. There have been no prior theoretical studies on systems where both 
reversal and rotation occur simultaneously in a large AC magnetic field. 
Consequently, comparisons with prior studies were performed under two extreme 
conditions. In the first case, high viscosities were assumed. Because reversal dom- 
inates rotation under these conditions, the results were compared with those reported 
by Carrey et al. 14 . In the second case, a high anisotropic field was assumed. Because 
rotation dominates reversal in this situation, the results were compared with the 
numerical simulations of Yoshida et al, where nonlinear Brownian rotational 
relaxation of magnetic fluids with a large excitation field was studied using the 
Fokker-Planck equation 13 . The results obtained from our simulation of reversal and 
rotation were consistent with those of earlier studies (Supplementary Figs. S1-S4). 
Therefore, we can now take a first step toward understanding the roles that rotation of 
a nanoparticle and reversal of its magnetic moment play together in large AC mag- 
netic fields. 
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